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Abstract 

The phonon dispersion curves of cubic BaTiOs have been computed within a 
first-principles approach and the results compared to the experimental data. 
The curves obtained are very similar to those reported for KNbC>3 by Yu and 
Krakauer [Phys. Rev. Lett. 74, 4067 (1995)]. They reveal that correlated 
atomic displacements along <100> chains are at the origin of the ferroelectric 
instability. A simplified model illustrates that spontaneous collective displace- 
ments will occur when a dozen of aligned atoms are coupled. The longitudinal 
interatomic force constant between nearest neighbour Ti and O atoms is rel- 
atively weak in comparison to that between Ti atoms in adjacent cells. The 
small coupling between Ti and O displacements seems however necessary to 
reproduce a ferroelectric instability. 

Keywords: BaTiOs, phonon dispersions curves, ferroelectric instability, den- 
sity functional theory, linear response. 
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I. INTRODUCTION 



Barium titanate (BaTiOs) is well known to exhibit three ferroelectric phase transitions: 
stable at high temperature in a perovskite cubic phase, its structure becomes successively 
tetragonal, orthorhombic and finally rhombohedral as the temperature goes down. Many 
works investigated and discussed the origin of its ferroelectric phase transitions |J. Amongst 
them, the most gratifying contribution was probably due to Cochran who associated the 
ferroelectric instability with the softening of a transverse optic phonon, and emphazised the 
connection between the structural instability and the lattice dynamics. 

Consequently to Cochran's work, a large amount of experiments have been performed in 
order to confirm the existence of a soft ferroelectric mode in BaTiOs. They include infra- 
red HQ] and Raman || measurements of the T phonon modes as well as various neutron 
diffraction data P-fLlfl. These experiments focused on the temperature behaviour of the 
soft phonon and were mainly concerned by the low frequency modes: only a few results are 
available for the higher energy vibrations. 

Simultaneously, theoretical phonon dispersion curves were deduced from fits to the exper- 
imental data using different shell models. Let us mention the pseudo-ionic model developped 
by Gnininvi and Bouillot [12[] or the rigid-shell model used by Jannot et al. [[LI]. These mod- 



els were however not particularly suited to describe ABO3 crystals, like BaTiOs. Migoni, 
Bilz and Bauerle |H| pointed out that the behaviour of the ferroelectric soft mode in the 
oxidic perovskites originates from the unusual anisotropic polarizability of the oxygen that, 
in turn, may be connected to hybridization between 0-2p and B-d states. Consequently, a 
more sophisticated "polarizability model" was introduced in order to include the spe- 



cific physical features of ABO3 compounds. The application of this model to BaTiOs was 



reported by Khatib et al. | 15fl . In their work, they obtained a full phonon band structure 



and investigated the temperature behaviour of the ferroelectric soft mode. However, their 
interesting results still remained at a semi-empirical level. 

Since a few years, theoretical advances have enabled one to determine the phonon dis- 
persion curves of solids within a truly first-principle approach. Recently, ab initio results 
have been reported for KNbOs [0 and SrTiOs []T7| . Similarly, in this paper, we investigate 
the lattice dynamics of BaTiOs. 

Our phonon frequencies compare well with the experimental data available. Moreover, 
our dispersion curves are very similar to those reported for KNbOs ||16|| , a perovskite material 
presenting the same sequence of phase transition as BaTiOs. in particular, the wave vector of 
the unstable phonon modes remains located in three (110) slab region of the Brillouin zone, 
emphasizing a chain-structure instability in real space fl6fl . This behaviour is illustrated 
with a simplified model. The range and anisotropy of the interatomic force constants is also 
discussed. 



II. TECHNICALITIES 



Our calculations have been performed in the general framework of the density functional 
formalism [18]. The exchange-correlation energy functional was evaluated within the local 



density approximation, using a polynomial parametrization of Ceperley-Alder [[HJ homoge- 
neous electron gas data. 
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We used a plane-wave pseudo-potential approach. The all-electron potentials were re- 
placed by the same extended norm conserving highly transferable pseudopotentials as in 
Ref. [^DJ. We have considered 5s, 5p and 6s as valence states to build the Ba pseudopoten- 
tial, 3s, 3p, 3d and 4s valence states for the Ti pseudopotential and 2s and 2p valence states 
for the O pseudopotential. The electronic wave function was expanded in plane waves up to 
a kinetic energy cutoff of 45 Hartree (about 6200 plane waves). Integrals over the Brillouin 
zone were replaced by a sum on a mesh of 6X6X6 special fc-points [2~I|. 

The dynamical matrix, Born effective charges and dielectric tensor were computed within 
a variational formulation of the density functional perturbation theory [ 2"3fl . First, 



calculations were carried out to determine the dynamical matrix on different meshes of q- 
points. Then, an interpolation was performed following the scheme proposed in Ref. |24|,|25|]. 
This technique takes properly into account the long-range behaviour of the dipole-dipole 
interaction which is separated from the remaining short-range forces owing to the knowledge 
of the Born effective charges and the optical dielectric tensor. 

An insight into the convergence reached on the phonon dispersion curves is reported in 
Fig. [I]. The frequencies deduced from the dynamical matrix at q — (.125, .125, .125) and 
q = (.375, .375, .375) are compared to those extrapolated from two different meshes of q- 
points: the first mesh (Ml) includes T (.0, .0, .0), X (.5, .0, .0), M (.5, .5, .0) and R (.5, .5, 
.5) points; the second mesh (M2) is the cubic mesh Ml to which the A (.25, .25, .25) point 
was added. We obtain a very good convergence with the M2 mesh. 



III. THE PHONON DISPERSION CURVES 

Our calculations are performed at the experimental lattice parameter of 4.00 A. This 
choice facilitates the comparison with the experimental data. Some indications on the vol- 
ume dependence of the phonon frequencies can be found in Ref. ||26|| , where the frequencies 



of the T phonons at different lattice constants have been compared. 

First, we present in Table | the Born effective charges and the optical dielectric con- 



stant fl27| , important ingredients of the present study since their knowledge allows to identify 
the long-range part of the interatomic force constants and makes the interpolation of phonon 
frequencies tractable. For Ba and Ti atoms, the effective charge tensor is isotropic. For O, 
the two independent components of the tensor correspond respectively to a displacement of 
the atom parallel (Zq,,) and perpendicular (Zq ± ) to the Ti-0 bond. (+7.32) and Zq^ 
(—5.78) are anomalously large [F2"U| , p8"| -|3"0]j with respect to the nominal ionic charges (+4 for 
Ti and —2 for O). This surprising phenomenon was explained recently in connection within 



dynamic changes of hybridization between 0-2p and Ti-3d orbitals PU| , |3"T| | . This specific 



feature is at the origin of a large destabilizing dipole-dipole interaction, connected to the 
anomalous mode effective charge [^,^] and the instability of the ferroelectric mode [26[ . 



Our computed optical dielectric constant (6.75) largely overestimates the experimental 
value (5.40) [|3l| , as usual within the local density approximation. This discrepancy will 
essentially affect the position of the highest longitudinal optic mode: when replacing the 
theoretical dielectric constant by the experimental value, its frequency at the V point changes 
from 631 to 696 cm -1 . At the opposite, the frequencies of the two other longitudinal modes 
at the T point are affected by less than 2 cm -1 . 
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The calculated phonon dispersion curves are plotted along high symmetry directions in 
Fig. ||. The T-X, T-M and T-R lines are along the <100>, <110> and <111> directions, 
respectively. The unstable modes associated to a negative curvature of the energy hypersur- 
face have imaginary phonon frequencies. The frequencies at the high symmetry points are 
reported in Table [TI[ 

Our result can be compared to the experimental data [|3|j6^p~Tj] . However, a complication 



arises from the fact that all the experimentally observed vibrational excitations have a real 
frequency while the computed unstable modes are obtained with an imaginary frequency. 
As the soft mode can be clearly identified by its symmetry, the associated experimental 
frequencies were removed from the comparison, for clarity. In the low frequency region, 
the presence of this additional soft mode may have slightly modified the frequency of the 
other modes. In spite of these difficulties we observe a good correspondence between our 
theoretical frequencies and the experimental data, specially for the acoustic modes for which 
a large variety of data are available. 

The ferroelectric phase transitions are driven by the unstable phonon modes. We are 
therefore mainly concerned by the analysis of these specific phonons. Two transverse optic 
modes are unstable at the T point: they correspond to a displacement of the Ti atom against 
the oxygen cage. The associated displacement eigenvector is equal to [<5(Ba) = —0.002, 
5(Ti) = -0.096, (J(Oi) = +0.158, 5(0 2 ) = 5(0 3 ) = +0.071] |3§. These two modes remain 
unstable all along the T-X line, with very little dispersion. One of them stabilizes along the 
T-M and X-M lines. Examination of the eigenvectors reveals that the unstable mode at the 
M (.5, .5, .0) point is polarized along the ^-direction: its displacement eigenvector is equal to 
[<5(Ti z ) = —0.130, 5(Oi )Z ) = +0.106] p4] . Both of the unstable modes become stable when 
deviating from the three T-X-M planes to the R-point. 

These features were also observed for KNbOa [Hj] and point out a marked 2D character 



of the instability in the Brillouin zone. This behaviour is more easily visualized in Fig. y 
where we show the frequency isosurface of the lowest unstable phonon branch corresponding 
to uj = 0. The region of instability, u; 2 (q) < 0, lies between three pairs of flat surfaces, 
that are parallel to the faces of the Brillouin zone cube. In other words, the unstable modes 
are contained in three perpendicular interpenetrating slab-like regions of finite thickness 
containing the T point. 



As highlighted by Yu and Krakauer [I6|, this behaviour corresponds to chain instabilities 
in real space. At the M-point, we have seen that there is a single unstable mode polarized 
along the z-axis and dominated by the Ti z and 01 z displacements. At this wave vector 
{iz — 0), the Ti and Oi atoms will be coherently displaced all along an infinite <001> 
chain. Going now from M to the R-point, the coherency of the displacement will gradually 
disappear and a finite length of correlation will be reached for which the phonon becomes 
stable. The finite thickness of the slab region of instability therefore corresponds to a 
minimum correlation length of the displacement required to observe an unstable phonon 
mode. From Fig. |3|, the length of the shortest unstable chain can be estimated to 4 a ce u = 
16 A |[j5|| . Note finally, the small dispersion of the unstable mode in the T-X-M plane 
suggests a small correlation of the displacements between the different Ti-0 chains. 
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IV. THE INTERATOMIC FORCE CONSTANTS 



In cubic BaTi03, we will see that the single displacement of a particular atom never 
leads to an instability: When one atom is displaced, a force is induced and brings it back 
in its initial position (see Tables |T| - [V|: the self-force on Ti and O is positive). However, 
its atomic displacement will simultaneously induce forces on the other atoms. It is only 
the additional displacement of some other atoms in this force field that can lower the total 
energy and produce an instability. The amplitude and the range of the interatomic force 
constants (IFC) associated to this mechanism can be analysed in order to clarify the chain 
instability pointed out in the previous Section. Moreover, the specific role of the dipole-dipole 
interaction (DD) can be separated from that of the short-range forces (SR) following the 
scheme proposed in Ref. J25|| . Our conventions on the interatomic force constants C ai p(K, k!) 
are such that the force F a (K) induced on atom k by the displacement At^k') of atom k' is 
given by: F tt («) = -C aj p(n, k') . At^k'). 

Let us first investigate the IFC with respect to a reference Ti atom along a Ti-0 chain 
(Table |ITT|) . As previously mentionned, we note that the self-force on the Ti atom is large 
and positive (+0.15215 Ha/Bohr 2 ). We observe also that the longitudinal IFC with the 
first neighbour O atom is surprisingly small (+0.00937 Ha/Bohr 2 ); moreover, it is positive. 
The analysis of the DD and SR contributions points out that these characteristics are the 
result of a destabilizing DD interaction, sufficiently large to compensate the SR forces. It is 
this close compensation which allows the displacement of Ti against the O atoms. Another 
insight on this balance of forces was already reported in Ref. p6| . Consequently to the very 
small total IFC, the Ti and O displacements might be relatively decoupled. 

At the opposite, the DD forces induced on the next Ti atom are negative: they will 
combine with the SR forces in order to produce sizeable coupling (—0.06721 Ha/Bohr 2 ). 
This mechanism is at the origin of the chain correlation of the Ti atomic displacements. By 
contrast, the transverse force on the first Ti neighbour is very small and confirms the small 
correlation of the displacements from chain to chain. 

The decay of the Ti— Ti and 0-0 longitudinal IFC with the interatomic distance can also 
be investigated. The results are reported in Table [TV] and [VJ. It is seen that the longitudinal 
IFC are anisotropic: they propagate essentially along the Ti-0 chain. This appears clearly 
for the SR part. For O, the DD contribution is also highly anisotropic due to the anisotropy 
of the Born effective charges. The anisotropy of the IFC is inherent to the chain correlation 
previously mentionned in this paper. 

V. THE CHAIN-STRUCTURE INSTABILITY 

From the knowledge of the IFCs previously reported, we can also investigate the energy 
surface of BaTi03 and illustrate the chain correlation highlighted in Section III. Let us 
consider that we have a bulk cubic crystal with the atoms frozen at their equilibrium position 
t ko . Then, we allow displacements of Ti and O atoms belonging to a [100] single Ti-0 chain 
of finite but increasing size. The total energy of this system will be given by: 

E(t k ) = E(t K0 ) + C 1j1 (k, k') At k At k , 
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where C is the interatomic force constant matrix and the sum on k and k' is restricted to 
the Ti and O atoms that are allowed to move. With the help of this equation, we can track 
the appearance of an instability in terms of the length of the chain of displaced atoms. An 
instability will correspond to a specific displacement pattern that lowers the total energy 
of the system: it will be associated to a negative eigenvalue of the restricted force constant 
matrix. 

In Fig. [|, we report the evolution of the lowest eigenvalue of the force constant matrix 
with respect to the length of the chain of moving atoms. Displacing only a single atom, the 
force induced on the Ti is larger than that on the O atom. With 3 atoms, we observe, at the 
opposite, that the Ti-terminated chain (Ti-O-Ti) is more stable than the O-terminated one 
(O-Ti-O): it points out the important role of the Ti— Ti interaction. The difference between 
Ti and O terminated chains will disappear progressively with the chain length. It is seen 
that an instability takes place for a chain longer than 10 atoms (5 unit cells). This is in 
close agreement with the correlation length estimated in the previous Section. It suggests 
that the behaviour of BaTiC>3 is already well reproduced when considering an isolated Ti-0 
chain of displacements. It confirms that the correlation between the different chains may 
play a minor role. 

Going further and frozing all the O atoms in such a way that only the Ti atoms are allowed 
to move along the chain, we can repeat the previous investigations. For this case, however, 
we do not observe any instability even for an infinite chain of correlated Ti displacements. 
This result aims to prove that the relatively weak coupling between Ti and O displacements 
still remains an important feature in the appearance of the structural instability. 

VI. CONCLUSIONS 

In summary, we have reported first-principles linear response calculations of the phonon 
dispersion curves of cubic BaTiC>3. Our results are in good agreement with the experimental 
data. Moreover, they are very similar to those reported for KNbO"3, a cubic perovskite 
material presenting a sequence of phase transitions analogous to BaTiC^. In particular it 
was seen that a ferroelectric instability takes place when correlated atomic displacements 
are allowed along <100> chains of finite length. This has been investigated with a simplified 
model and the length of the shortest unstable chain estimated to a dozen of atoms. The 
interatomic force constants are anisotropic and propagate essentially along the Ti-0 chains. 
Their analysis has emphasized the correlation of the Ti displacements. It has shown that 
the Ti and O atomic displacements are only weakly coupled. This small coupling remains 
however an important ingredient to reproduce the ferroelectric instability. 
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FIGURES 



FIG. 1. Convergence achieved on the calculated phonon dispersion curves of cubic BaTi03 
along the T-R line. The open symbols correspond to g-points included in the Ml (circle) and M2 
(circle+square) meshes used to extrapolate the curves (Ml: dotted lines; M2: full lines). The 
filled symbols are associated to points not included in the mesh: they illustrate that a satisfactory 
convergence is obtained with the M2 mesh. 

FIG. 2. Calculated phonon dispersion curves of cubic BaTiOs at the experimental lattice con- 
stant. The theoretical result shows a reasonable agreement with the experimental data: (•) Ref. 
[3], (o) Ref. [6], (+) Ref. [7], (□) Ref. [8], (x) Ref. [9], (V) Ref. [10], (A) Ref. [11]. 

FIG. 3. Zero-frequency isosurface of the lowest unstable phonon branch over the Brillouin zone. 
T is located at the center of the cube. The mode is unstable in the region between the nearly flat 
surfaces. 

FIG. 4. Lowest eigenvalue of the restricted force constant matrix associated to atomic displace- 
ments along a finite Ti-0 chain of increasing size. 
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TABLES 



TABLE I. Born effective charges and optical dielectric constant of cubic BaTi0 3 at the exper- 
imental volume [27]. 



Experiment Present Zhong et al. c 







+2.9° 


+2.74 


+2.75 






+6.7 a 


+7.32 


+7.16 






-2.4 a 


-2.14 


-2.11 






-4.8 a 


-5.78 


-5.69 


f oo 




5.40 6 


6.75 




a Reference 


28 








b Reference 


32 








c Reference 


29 









TABLE II. Computed phonon frequencies (cm 1 )of cubic BaTi03 at T, X, M and R. Symmetry 
labels follow the convention of Ref . [33] . 



q-point 


label 


frequency 


label 


frequency 


r 


Li 5 (TO) 


219 i 


r 25 


281 




Lis (A) 





Ti 5 (LO) 


445 




Tis (LO) 


159 


Li 5 (TO) 


453 




Li 5 (TO) 


166 


Ti 5 (LO) 


631 


X 


x 5 


189 % 


x 3 


322 




X 5 / 


104 


X 5 / 


330 




X2' 


146 


x 5 


421 




x 5 


194 


Xi 


517 




Xi 


260 


X2' 


627 


M 


M 3 , 


167 i 


M 5 


344 




M 2 / 


103 


M 2 


354 




My 


104 


My 


435 




M 3 


208 


Mi 


456 




Mb/ 


270 


M 4 


683 




M 3 , 


333 






R 


Rl5 


128 


R25' 


386 




R25 


182 


Rl5 


414 




Rl2' 


314 


R2' 


717 
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TABLE III. Longitudinal (||) and transverse (_L) interatomic force constants (Ha/Bohr 2 ) with 
respect to a reference Ti atom (Ti(0)) along the Ti-0 chain of cubic BaTiC^. 



Atom 


Total force 




DD force 


SR force 


Ti{0) 


+0.15215 




—0.27543 


+0.42758 




+0.00937 




+0.23247 


-0.22310 


Ti\\(2) 


-0.06721 




-0.03680 


-0.03041 


0||(3) 


+0.01560 




+0.00861 


+0.00699 


Ti|| (4) 


-0.00589 




-0.00460 


-0.00129 


o±(i) 


-0.02114 




-0.04298 


+0.02184 


Ti±(2) 


+0.00751 




+0.01840 


-0.01089 


TABLE IV. Ti-Ti longitudinal interatomic force constants (Ha/Bohr 2 ) with respect to a refer- 


ence Ti atom at (.5, .5, 


.5). 








coordinate 


distance 


IFC 


DD part 


SR part 


(.5, .5, .5) 


0.0000 


+0.15215 


-0.27543 


+0.42758 


(-.5, .5, .5) 


7.5589 


-0.06721 


-0.03680 


-U.UOU41 


(-.5, -.5, .5) 


10.6899 


-0.01114 


-0.01301 


+0.00187 


(-.5, -.5, -.5) 


13.0924 


-0.00643 


-0.00780 


+0.00065 


(-1.5, .5, .5) 


15.1178 


-0.00589 


-0.00460 


-0.00129 


TABLE V. 0-0 longitudinal interatomic force constants (Ha/Bohr 2 ) with respect to 


a reference 


atom at (.5, .5, .0). 










coordinate 


distance 


IFC 


DD part 


SR part 


(.5, .5, .0) 


0.0000 


+0.12741 


-0.35322 


+0.48062 


(.5, .0, .5) 


5.3450 


-0.02838 


-0.03367 


+0.00529 


(-.5, .5, .0) 


7.5589 


-0.00190 


-0.00314 


+0.00124 


(.5, .5, -1.0) 


7.5589 


-0.03212 


-0.02295 


-0.00918 


(-.5, .0, .5) 


9.2577 


-0.00183 


-0.00289 


+0.00106 


(-.5,-.5, .0) 
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